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We propose an experimental procedure to directly measure the state of an electromagnetic field 
inside a resonator, corresponding to a superconducting transmission line, coupled to a Cooper-pair 
box (CPB). The measurement protocol is based on the use of a dispersive interaction between the 
field and the CPB, and the coupling to an external classical field that is tuned to resonance with 
either the field or the CPB. We present a numerical simulation that demonstrates the feasibility of 
this protocol, which is within reach of present technology 
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Superconducting electrical devices have been experi- 
mentally proven to be serious candidates for the realiza- 
tion of quantum information processing tasks . Coher- 
ent control and near unit-visibility Rabi-oscillations 2, 
3], coupling of two CPB-qubits jj] and the implemen- 
tation of conditional gates y| are striking experiments 
that demonstrate the high level of control achieved on 
these systems. Furthermore, a scalable architecture for 
quantum computation has already been envisioned @. 

On the other hand, recent demonstrations of Jaynes- 
Cummings-like dynamics between a CPB-qubit and the 
quantized mode of a superconducting transmission line 
resonator (which acts as a quasi-lD cavity) [3, have 
shown that many of the tools originally developed within 
the context of quantum optics can now be extended to 
solid state physics. Once coherent control and complete 
characterization of quantum states have been achieved 
at the qubit level, it is natural to attempt such levels 
of control for the electromagnetic field generated by the 
transmission line. For its characterization one could, in 
principle, make use of the well-known homodyne and 
heterodyne detection techniques. But, since the field 
we would like to characterize is inside a resonator and 
consists of a few photons, implementation of those tech- 
niques turns out to be a non-trivial task. Homodyne 
detection has been proposed for characterizing the state 
of the field leaking out from a tridimensional cavity in § , 
and for a one-photon field leaking out of a 1-D cavity in 
0- Nevertheless, it is very difficult to apply this proce- 
dure to high-finesse cavities containing weak fields, since 
one would have to distinguish a still weaker leaking field 
from the noise in the detector. Furthermore, unavoidable 
absorption losses may lead to poor reconstruction of the 
state of the intracavity field, as pointed out in ^(| • 

To overcome these issues we propose here an experi- 
ment to directly measure the Wigner function of the 
electromagnetic field inside a superconducting transmis- 
sion line resonator coupled to a CPB-qubit, via the mea- 
surement of the latter's populations. The Wigner func- 
tion contains all the information about the state of the 



field, and is a useful tool for studying the decoherence- 
induced quantum-to-classical transition, as it provides us 
with a phase-space representation that can be compared 
to classical probability distributions 01 ■ For a single 
mode of the electromagnetic field, it is defined in terms 
of the respective density operator p as [l3T |: 

W(a) = (l/ 7 r)Tr[ / 5i)(a)PLi~ 1 (a)]. (1) 

Here, D(a) = e aat ~ Q a is the field displacement oper- 
ator, which takes any coherent state \[3) to |/3 + a), up 
to a phase factor, and P = e~ l7r ° ta is the parity opera- 
tor, which multiplies a Fock state \n) by a factor (— l) n ; 
a and are respectively the photon annihilation and 
creation operators of the mode. The displacement oper- 
ator can be operationally implemented, in a cavity QED 
(cQED) setup ^4|, by injecting a coherent field with com- 
plex amplitude a into the cavity. A protocol for the di- 
rect measurement of the Wigner function was first pro- 
posed in [l5| and later experimentally carried out in [HI 
for the microwave field inside a 3-D high-quality-factor 
(Q) cavity. It involves injecting a microwave field (com- 
plex amplitude a) into the cavity, so as to displace the 
field to be measured, and then sending an atom with 
two of its levels, |e) and \g), interacting dispersively with 
the displaced field. The atom is prepared in the state 
(|e) + |<?))/V2, and, after leaving the cavity, is submitted 
to a classical field, so that its state undergoes a ir/2 rota- 
tion [je) -> (|e) + \g))/2, \g) -> (-|e) + \g))/2]. Then the 
atomic population is measured. The difference between 
the probabilities of finding the atom in states |e) and \g) 
is proportional to the value of the Wigner function of the 
cavity field at the point —a in phase space. 

It is not possible however to apply this protocol to the 
system here considered, since in this case the atom (CPB- 
qubit) is always inside the cavity and its interaction with 
the field cannot be turned off. Nevertheless, we show here 
that it is still possible to directly measure the Wigner 
function of the electromagnetic field in a superconducting 
transmission line, via the Copper-pair box qubit. Our 
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method could also be applied to other systems involving 
the interaction of a qubit with a resonator 

The system under consideration consists of a meso- 
scopic superconducting island (see Fig. capacitively 
coupled to the quantized field mode of a transmission line 
of length L (see Fig. ^>). Details of this system can be 
found in Q. The CPB Hamiltonian is given by 0: 



H CPB = 4E c (n - n g f - Ej cos 9 . 



(2) 



where h is the number operator corresponding to the 
Cooper-pair charges in excess on the island, and is the 
average phase drop along the junctions (n = —id/d(Q)). 
Each junction is characterized by a capacitance Cj and 
a Josephson energy Ej Q . The effective Josephson cou- 
pling Ej = 2Ej cos(7T$/$ ) can be changed through 
an applied static magnetic flux <J> ($ = h/2e is the 
flux quantum and e the electron charge). The charg- 
ing energy is E c = e 2 /2Cs (Ce = C' g + 2Cj) and the 
gate charge is n g = C g V g /2e, which can be tuned by 
the dc part of the potential gate V g . The coupling to 
the quantized field mode of the transmission line, of fre- 
quency u>, is taken into account through the quantum 
part of the gate voltage, i.e. V g = V g c + v, where v = 

y^huj/Lc (aJ + a) (c is the transmission line capacitance 
per unit length), with a being the annihilation opera- 
tor for the transmission-line mode. Finally, the Hamilto- 
nian of the CPB-resonator system is obtained by adding, 
to Eq. (J2J, the Hamiltonian of the oscillator mode, i.e. 
H os = Ticua^d. In the charge regime, i.e Aqap ^S> Eq ^S> 
AEj (Aqap is the superconductor gap) the CPB can 
be treated as an effective two-level system |lj of tran- 
sition frequency uj = ^Ej + [AE c (l - 2n g ic )] 2 /h, with 

rig C = C g V g c /2e. Within this regime, when n g c is around 
1/2 and in the rotating-wave approximation, the Hamil- 
tonian of the composite system reduces to that of the 
Jaynes-Cummings (JC) model with coupling Rabi fre- 
quency g = {eC g / C-z) \J Tito / Lc. In order to drive the 
composite system the transmission line is coupled ca- 
pacitively (C Q ) to an external classical microwave field, 
of frequency u>d and slowly varying complex amplitude 
e(t) = en{t)+iei{t), whose effect can be modeled through 
the driving Hamiltonian Hd = h[e(t)d^ e~ luJdt + H.c.]. A 
second-order perturbative calculation, in the dispersive 
regime \g\Jn + 1/A| -C 1 (n is the mean photon number 
and A = luq — lu is the detuning between the cavity and 
the CPB-qubit) yields, for the total system dynamics (in- 
cluding the driving), in a reference frame rotating with 
the driving field frequency, the effective Hamiltonian 

Hlf f = n(u-u d )tfa+h [e(t)a f + H.c,]+(h/2)n-& , (3) 
where 



fl = pg/A)e R (t),-(2g/A)6i(t),LJo - Lu d 
+ (g 2 /A)(2^a+l)} , 
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FIG. 1: a) The Cooper-pair box: a superconducting island with n 
Cooper pairs in excess (dashed area) with two Josephson Junctions 
(crossed boxes), b) Resonator-CPB-qubit composite system: the 
central superconducting waveguide (black) is coupled capacitively 
(Co) to other two superconducting lines (dashed), through which 
the classical driving fields are pumped. 



and a = (<j x , (t v ,(t z ) are the Pauli matrices. 

Hamiltonian © generates a displacement of the field 
but also induces a rotation on the qubit Bloch vector, the 
simultaneity of both operations coming from the impos- 
sibility of switching off the cavity-qubit interaction. 

The measuring protocol consists of first encoding in- 
formation contained in the initial field state, p, into the 
CPB-qubit populations, which are then measured. The 
encoding process is divided into four evolution steps: 1) 
a coherent displacement of the field; 2) a 7r/2-pulse on 
the CPB-qubit; 3) a dispersive evolution without driv- 
ing; and finally 4) another 7r/2-pulsc on the CPB-qubit. 
The displacement of the field as well as the pulses on the 
qubit are driven by the external classical microwave field, 
whose complex amplitude e(t) — e and frequency tod are 
the parameters under control. 

In the first step we set ujd = u> — <? 2 /A and e = er>, 
with |erj| g + A 2 /2g. This last condition implies, from 
@ and (@J, that we can neglect the components of the 
rotation axis Q perpendicular to the z-direction. So, if 
initially the qubit is in its ground state |<?), it remains 
there, and the evolution operator in a reference frame 
rotating with frequency u — g 2 /A is given by Ui(tD, 0) = 
D(a), where a = —ieotry, and tr> is the pulse duration. 

In the second step we set the driving frequency cud = 
and e = |e 7r / 2 |e 1 ^ 1 , during a time t v / 2 = nA/ Ag\e^/ 2 \- 
Then, according to l@J, each Fock component of the 
state suffers a 7r/2-rotation of the Bloch vector about 
an axis whose z-component, {g 2 / A)(n + 1/2), depends 
on the photon number n. This component can be ne- 
glected under the condition \e„/ 2 \ ^> g(n + 1/2) n <C 
\ e K/2\/g ~ 1/2, where n is the mean number of photons 
in the state [l8|- In this case, the evolution operator is 

given by U 2 (t D + t n/2 ,t D ) = e lAatat '/^K /2 )i? s (V2) 
in the representation rotating with frequency loq (which 
will be used from now on) . This essentially consists of a 
displacement, D{a^/ 2 )i 01 the field state by an amplitude 

iAt D r t "/2+*-D , , „-iAt, 



' "~ tD /to™ e V2 e ~ iAtc ^ and the rotation, 



(4) /?fj (7r/2), of the qubit state by an angle it/2 about an 
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axis 17 in the equatorial plane of the Bloch sphere. 

In the third step we switch off the driving field, e = 0, 
and let the system evolve freely during a time tp = 
nA/2g 2 according to U 3 (t D + t„/ 2 + t P ,t D + t^j-x) = 

■ TV A 2 ~ \ ~ t 1 

e 2 s' 2 a a e~ l ^^ a a +2) a ' . Because the qubit is already in a 
superposition of the upper |e) and lower \g) states, this 
is the step where the field and the qubit get entangled, 
which is crucial for the transfer of information from the 
field to the qubit. The fourth and last step in the en- 
coding process is another 7r/2-pulse on the qubit, with 
e = |e Tr /2|e i ^ 2 , corresponding to an evolution operator 

U 4 analogous to C/ 2 - 

Collecting all the steps of the encoding process we 
get the total evolution Ut = U4U3U2UU1, where U = 

switches from a frame rotating 



i(A+g 2 /A)(a t a+o' i /2)t D 



g /A to one rotating with frequency 



with frequency u> 
uj . Thus, (a z ) 
straightforwardly calculated, yielding 

-fa)pD(^e 



P e - P g = TT[a z U T \g)(g\pU±} is 



p 



P g = Tr sinf 



D(0) , (5) 



where /3 = a + (2|e 7r/2 |/A) sin(t 7r/2 A/2)e~^, with a = 
-ie D t D and tj> = (t D + t 1l/2 /2)A + g 2 t D /A + vr/2 - fa. 
When fa — fa — tt/2 we see from Eq. (QJ that (<J Z )/^ 
yields the value of the initial-field Wigner function at 
the point —j3. If the first 7r/2-pulse is chosen so that 
£„./ 2 A/2 = ran, m integer, then (3 = a. By repeating the 
experiment for different t]js one can scan the whole phase 
space and thus fully reconstruct the quantum state. 

Table |I| displays a comparison between the parame- 
ters reported in Q and the optimal parameters that we 
propose. The performance of the protocol for these two 
sets of parameters were tested by a numerical simula- 
tion, where each step of the protocol was carried out 
evolving the system with the exact JC-model plus the 
driving Hamiltonian. The Wigner function of the field 
state was obtained using Eq. (J5J where the probabili- 
ties P e and P g were calculated for the final entangled 
state of the system (see Fig. 0). With the first set of 
parameters, the whole measurement protocol takes ap- 
proximately 100 ns, which is less than the cavity lifetime 
k _1 and much smaller than the atom lifetime 7 -1 . More- 
over, the condition \efa -C g + A 2 /2g w A 2 /2g = uj, 
used in the first step of the protocol, is well satisfied. 
On the other hand, we see that the condition for the 
7r/2 rotation is not comfortably met, since it is not true 
that n <C |&7r/2|/5 — 1/2 = 4.5 for the states of inter- 
est. In fact, because the mean photon number increases 
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TABLE I: Experimental parameters. A, g, cd, and e n /2 are 
expressed in units of the transmission line frequency. The first 
set of parameters is reported in ; the second is the proposed 
one. K _1 is the cavity lifetime and 7 -1 is the atom lifetime. 
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FIG. 2: Wigner function W(a) on the real axis a = x for different 
states of the cavity field: (al) and (a2) Vacuum state; (bl) and (t>2) 
Fock state with n = 1; (cl) and (c2) Schrodinger-cat state oc (|ao) — 
— cto))i with 00 = 2. The full line is the theoretical value and 
the crosses stand for the values obtained by a numerical simulation 
of the measurement protocol (see text for details). The numerical 
simulations on the left-hand (right-hand) side were performed with 
the first (second) row of parameters in Table The transmission 
line frequency is taken as w = 2n X 10 GHz. 



after the displacement of the field in the first step, i.e. 
n — * n + 25Re(a*(a)) + |a| 2 , the conditions for the 7r/2 
rotation and the dispersive regime approximation are vio- 
lated for greater values of |a| (see the numerical points on 
the tails of the graphs on the left in Fig.J2J)- During the 
7r/2 rotation of the CPB-qubit the field is also displaced 
to (3 = a + (2|e OT/2 |/A)sin(tA/2)e-^, so the displace- 
ment |/3 — a I attains its maximum value at the middle 
of the rotation, violating the dispersive-regime approxi- 
mation for the first set of parameters. These considera- 
tions imply that the accuracy of the method is worse for 
the tails of the Wigner functions, since probing them re- 
quires larger displacements of the cavity field. The poor 
accuracy in Fig. 2(cl) is due to the contribution of high- 
n Fock states. Finally, one should consider that with 
= 160 ns decoherencc effects may become apprecia- 
ble at 100 ns (duration of the measurement protocol in 
this case) for fields with average photon number larger 
than one. Thus, a higher-Q cavity should be required. 

With the detuning proposed in the second set of pa- 
rameters the total duration of the protocol would be 
about 300 ns, for a cavity lifetime k _1 ~ 1000 ns. The 
difference of one order of magnitude between the val- 
ues for the damping times in the two sets of parame- 
ters in Table I is easily overcome with present technol- 
ogy. Better lifetimes could be achieved, for example, by 
decreasing the value of the capacitances Co- On the 
other hand, |e w / 2 | increases with the power of the ex- 
ternal radio-frequency source. The improvement in the 
reconstruction of the Wigner function for this second set 
of parameters is displayed in the graphs on the right- 
hand side of Fig. [3 for the vacuum state, the Fock state 
with n = 1, and a Schrodinger- cat state. The main 
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impact of using the new parameters is a great reduc- 
tion of the errors in the tails of the Wigner functions. 
Decoherence, not taken into account in our simulations, 
would further limit the maximum number of photons in 
the states characterized. Temperature effects are negligi- 
ble for typical experimental values (T = 100 mK in fj). 
Indeed, for T = 100 mK, and given that for our param- 
eters we have hui/k w 480 mK, the thermal occupation 
number is n < 0.008. The reconstruction of states with 
more photons would demand higher Q's, but the proto- 
col would remain the same. We note that the experiment 
could be used to continuously monitor the loss of coher- 
ence of the field. 

As for the measurement of the qubit population, a dis- 
persive quantum non-demolition scheme was carried out 
in 0. However, for the higher Q value considered here, 
this technique would take a time of the order of the qubit 
lifetime. The direct measurement of the qubit population 
could be accomplished in this case by coupling to the 
qubit a single electron transistor (SET) device, which is 
able to detect charge differences of the order of e, as de- 
scribed, for example, in 19]. The influence of the SET 
on the qubit dynamics can be minimized by turning the 
device on only at the moment of measurement, as dis- 
cussed in 0] • The presence of an extra superconducting 
lead connecting the SET to the qubit should not increase 
significantly the decoherent effects on the cavity field al- 
ready introduced by the presence of the Cooper-pair box. 

Experimental testing of this protocol would require the 
preparation of simple field states. Coherent states are 
prepared by displacing the initial state \g,0), which is 
accomplished by setting the driving parameters lu c i = uj — 
g 2 / A and e = eo, as in the first step of our protocol. This 
will take a time tr> = \a\/\eo\ to be carried out, where a 
is the amplitude of the coherent state. The generation of 
a Schrodinger cat-like state, i.e. (\ao) ± e lip \ — a ))/N, 
where N is a normalization factor, is contained implicitly 
in the protocol described in this paper, since after the 



whole evolution stage, and before the qubit population 
is measured, the final entangled state of the system is of 
the type (1/V2)(|a ) + | - a )) <8> |e) +(1/V2)(|a ) - 
e l(p | — ao)) ® Iff)- Measuring the qubit population would 
project this state onto a coherent superposition of two 
coherent states. The time required to generate this state 
is the same as for our protocol. 

Another example of interest is the one-photon Fock 
state: beginning with the system in the state |<?,0) one 
applies a 7r-pulse on the qubit, setting LOd = u>q and e = 
during a time = nA/2g\e 7r \. Choosing t n A/2 = mir, 
with m integer, the state |e, 0) is prepared. Then, we 
tune the qubit frequency luq into resonance with the cav- 
ity mode (by changing the magnetic flux $) and let the 
system complete a Rabi oscillation (tR a = ir/2g), so the 
resulting state is \g, 1). Next we change the magnetic flux 
again to take the system back to the dispersive regime. 
This method would require rapid switching of the flux 
(less than 1 ns), a challenge for present experiments. 
Alternatively, time-dependent magnetic fluxes could be 
used to tune the qubits into and out of resonance with the 
cavity field, as recently suggested in |2l| . For the second 
set of parameters in Tableland \e w \ = 0.3 u> (m = 15) the 
total time for this process would then be t^a + t-n ~ 10 
ns. 

In conclusion, we propose here an experiment to com- 
pletely characterize the electromagnetic field of a quasi 1- 
D superconducting transmission line resonator, which is 
always interacting with a CPB-qubit, by directly measur- 
ing its Wigner function. This is carried out by measuring 
the CPB-qubit population after the application of a se- 
ries of driving pulses, induced by an external microwave 
field coupled to the resonator. Our numerical simula- 
tions, for realistic parameters, show that this method is 
within reach of present experimental setups. 
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